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Abstract 

It is generally accepted that the observed CMBR dipole arises from the 
motion of the local group relative to the CMBR frame. An alternative inter- 
pretation is that the dipole results from an ultra-large scale (A > lOOc/i^o) 
isocurvature perturbation. Recently it was argued that this alternative pos- 
sibility is ruled out. We examine the growth of perturbations on scales larger 
than the Hubble radius and in view of this analysis, we show that the isocur- 
vature interpretation is still a viable explanation. If the dipole is due to 
peculiar motion then it should appear in observations of other background 
sources provided that they are distant enough. 

I. INTRODUCTION 

The dipole moment is the most prominent feature in the Cosmic Microwave Background 
Radiation (CMBR) anisotropy [1]. The dipole, which is larger by two orders of magnitude 
than all other multipoles, is generally accepted to result from the earth's motion relative to 
the CMBR frame. The main purpose of this paper is to emphasize the fact that the origin 
of the dipole is not necessarily a Doppler effect. This idea has already been put forward 
by a few authors ( [2], [3], [4]). However, it was recently claimed that these arguments 
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were wrong [5] . We show here that a large scale isocurvature model for the dipole is a viable 
alternative to the Doppler origin. Current observations of the CMBR dipole and quadrupole 
are consistent with this possibility. We examine other potential implications of this scenario 
and suggest observations that would confirm or rule out the Doppler origin of the dipole. 

Our starting point here will be the paper by Paczynski and Piran [3], hereafter denoted 
PP. This paper is based on a Tolman-Bondi model (spherical symmetry and dust), which 
contains a (gravitationally negligible) spherical distribution of radiation. It is shown that 
a non centered observer can measure a significant dipole due to such a radially varying 
specific entropy (i.e. the ratio of the number density of photons to the number density of 
baryons). These results were obtained by integrating numerically the light geodesies in the 
Tolman-Bondi model. 

The phenomenon described in PP can appear to the reader a bit artificial by the choice of 
a very particular space-time model and the results are not intuitive in view of the complicated 
numerical integration involved. Moreover, there was a recent claim [5] that the results of PP 
are wrong and that it is impossible to obtain a dipole far larger than the quadrupole from 
either isocurvature or adiabatic perturbations. It is, therefore, the purpose of this paper to 
demonstrate that the phenomenon described in PP is on one hand true and on the other 
hand is more general than it seems at first glance. To do so we show that the results of PP 
can be obtained in the context of a linearly perturbed Friedmann- Robertson- Walker (FRW) 
universe, within which the computation of the dipole and of the quadrupole can be carried 
out analytically. These results depend, in fact, only on one crucial argument: the presence 
of very large scale isocurvature perturbations (by very large scales, we mean scales far larger 
than the Hubble radius today). It is essential to stress that the modes that contribute to 
the observed CMBR dipole and quadrupole anisotropy are much larger than the Hubble 
radius at last scattering. Therefore, we will focus our analysis only on the evolution of the 
perturbation modes outside the Hubble radius. 
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We also wish to answer to another objection which could have been made against the 
model of PP, namely the fact that they consider our universe only in the phase of matter 
domination and that they add an ad hoc isocurvature perturbation at the time of the last 
scattering. A question of interest is whether some primordial isocurvature perturbations can 
survive during the evolution of the universe and be sufficiently important at the time of last 
scattering to produce effects comparable to those in PP. We show that this indeed the case 
by considering the influence of a pure isocurvature primordial perturbation on the dipole 
and quadrupole moments. 

The plan of this paper is the following. In the section 2, we introduce the concept of adi- 
abatic and isocurvature linear perturbations in a flat FRW background and we rederive the 
equations governing their evolution. In section 3, we give the expression for the anisotropy 
of the CMBR. In section 4, we make the connection between the Tolman-Bondi model used 
by PP and our cosmological perturbations approach. Finally we summarize in section 5 the 
observational implications of these results. 

II. ADIABATIC AND ISOCURVATURE PERTURBATIONS 

There are several formalisms for dealing with cosmological linear perturbations. The 
oldest is due to Lifschitz [6] and uses the synchronous gauge. Another is the so-called 
gauge-invariant formalism of Bardeen [7], which employs arbitrary gauge and constructs 
gauge invariant quantities out of linear combinations of the perturbations. Finally there 
is also a covariant approach of cosmological perturbations, pioneered by Hawking [8] and 
developed recently by several authors [9] (see also references in [10]). We use here this 
latter formalism which we find more convenient: it employs quantities with a clear physical 
meaning, and in particular it provides a direct definition of the peculiar velocity, which turns 
out to be useful in interpreting the Sachs- Wolfe effect (see [If]). 

We begin by reviewing this cosmological perturbation theory for a single fiuid. Consider 
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a space-time, endowed with a metric (/^j,, filled with a perfect fluid with a number density, n, 
an energy density, /?, a pressure, p and a four- velocity u^. One defines the comoving gauge 
as a particular foliation of space-time into hypersurfaces that are orthogonal to the matter 
fiow u''. This foliation is identified with the preferred foliation of the FRW spacetime (which 
we take to be fiat for simplicity) representing the homogeneous approximation of the real 
spacetime. Note that such a foliation always exists in the linear approximation (whereas, in 
general, it requires a vorticity free fiow). 
We define a local Hubble parameter by 

3H = Vy. (1) 

Using this definition, the local conservation of matter, 

V^(nu^) = 0, (2) 

can be written as 

^ = -?>H{p + p), (3) 

where dj dr is the derivative along the fiow lines u^V^. One can, then, write the Raychaud- 
huri equation, ignoring the terms involving the shear and the vorticity which are second 
order in the perturbations, as 



dr 6 6 p + p 



H'--^ip+m-o—^- (4) 



The operator, D^, stands for D^D^ where is the covariant derivative projected on the 
hypersurface orthogonal to the fiow u^, i.e. = [g'j^ -\- u^u'^)\^^. Equations 3 and 4 can be 
combined to a single equation governing the time evolution of 6 = Sp/ p: 

^6 + h(2 + 3cI- 6w) ^6-Ih'(i + Sw- 3w' - Gc^) 6 = (5) 



' V ' « J dt 2 ^ '> p 

where t is the time parameter of the comoving hypersurfaces. It is related to the proper 
time by ( [12]) 



The previous analysis deals only with a single perfect fluid. One can extend this treat- 
ment to several uncoupled perfect fluids (see [f3] and [f2]), each specified by its four- velocity 
M^, its energy density pa and its pressure Pa- We are interested only in first order deviations 
from the FRW configuration where all the fiuids have a common four- velocity. For each fiuid 
one can define a Hubble parameter Ha by an equation similar to (f). Moreover each fiuid 
satisfied a conservation equation similar to (3) since the fiuids are decoupled. It can then 
be shown [12] that the Raychaudhuri equation for each fiuid becomes, ignoring higher order 
perturbative terms, 

^^H'^^{p + ?>p) = -i-^ + ^'^"" pa. (7) 

dr " 3 3 Pa + Pa Pa + Pa 

By linearizing the conservation equation and the Raychaudhuri equation, one finds a system 
of coupled first order differential equations, 

d c nTT n ^P o /I , \ e TT , n TT ( 2 



and 



-8a = ?>HJa— -3(1+ Wa) SHa + SH^ " 6^, 

at Pa ^ ^ 



at 3 3 Pa + Pa Pa+ Pa 

where one has defined 

r. Pa ~\~ Pa Pa /i n\ 

0a = ■ , Wa = — . (10) 

P + P pa 

A dot denotes the time derivation for the homogeneous background quantities. 

If one considers only two fiuids, it is useful to introduce the perturbation in the ratio of 
the number density 



and the total density perturbation 

S = ^^1±^ = ^S, + ^S,. (12) 

Pl+ P2 P P 

A general perturbation can be described by the pair ((^1,(^2) or alternatively by the pair 
(S^S). An adiabatic perturbation satisfies S = and an isocurvature perturbation satisfies 
6 = 0. Unfortunately, these conditions are not invariant with time and a perturbation that 
begins as an isocurvature perturbation generates an adiabatic component and vice versa. 
This was the origin of the claim of [5] that the analysis of PP is wrong. However, as we 
show later, and [5] fail to realize, for perturbations larger than the horizon, if ~ initially 
S will remain very small with respect to 6 and vice versa. In particular this decomposition 
is meaningful for primordial fiuctuations. 

We now rewrite equations (8) and (9), adapted to the variables (61^62)^ as evolution 
equations for the quantities 6 and S. Using (8) and Wa = —'iH{c\ — ti'a)(l + Wa) (where 
c\ = Pal Pa) , one finds that 

^S = -?,{8H^-8H2). (13) 
at 

Differentiating this equation and using (9), we obtain a second order differential equation 
for S\ 

where = c\62 + c\6i. Finally the equation for the evolution of the density perturbations 
8 follows from (5) when one expresses 8P in terms of 8 and S: 

^ + h[2 + 3cI-6w)j^- '^H' [i + 8w- 3w' - Gc^) 8 = clD^8 + (c^ - cDOrO^il + w)D^S. 

(15) 

The physical situation, which we consider from now on, is the case where fiuid 1 is 
pressureless {c\ = 0) and fiuid 2 is radiation {c\ = 1/3). We define apq as the scale factor at 



the matter-radiation transition, i.e. when the energy densities of the two fluids are equaL 
Then, by using pi/p2 = a/aeq^ one has 

eJ = ifl + 5AV\ ,; = ifl4^y\ ,16) 



3 \ 4 ttpq I ^ 3 V 3 a 



au' 



d 



w = \(l + ^)\ {cl- cDOMl + w) = -\(l + l^)\l + ^y\ (17) 

Note that with this choice of fluids 1 and 2 our deflnition of S corresponds to the opposite 
of the variation of the speciflc entropy, i.e. the ratio of the number density of photons to 
the number density of dust (denoted S in PP). 

It is convenient to introduce the Fourier decomposition of the perturbations according 
to the deflnition 

S, = J (^^-"•"^(-)' (18) 

with a similar deflnition for 6^. The evolution equations for the Fourier modes are simply 
equations (14) and (15) modifled with the substitution of —k'^/a'^ in place of D^. Fach 
Fourier mode evolves independently of all the other modes, and one can thus study each 
mode individually. In standard cosmology the matter satisfles the strong energy condition 
3p + /O > and the comoving Hubble radius [aH)~^ increases with time as the universe 
expands (the inverse is true during inflation). This implies that any given Fourier mode 
was outside the Hubble radius at sufflciently early time. Consequently, it is traditional, 
in standard cosmology, to deflne the initial conditions for the perturbations during the 
radiation dominated era at a stage when the relevant mode was outside the Hubble radius, 
i.e. when k <^ aH . During this stage, the r.h.s. of equation (14) can be neglected and the 
corresponding solutions are ~ const, and ~ ln(t). The second solution is singular in 
the past and can be ignored. We see that primordial isocurvature perturbations are constant 
in time when they are outside the horizon: 
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(19) 



where the superscript 'p' denotes the primordial value. This remains true even after the 
transition between radiation domination and matter domination, as long as the modes are 
outside the Hubble radius (see Figure 1). 

The isocurvature mode generates an adiabatic perturbation. We have integrated nu- 
merically the evolution of the density mode 6k produced by a pure isocurvature primordial 
perturbation [6k = 0). We find that 6k grows. However, as long as the mode remains outside 
the Hubble radius, the value of 6k is small with respect to the corresponding value of Sk- 
During the radiation era (see [10]), 



The amplitude of the density perturbation mode continues to grow after the radiation-matter 
transition, but as long as the mode is outside the Hubble radius, the amplitude is bounded 
from above by the asymptotic limit (see Fig. 2 and also [10]) 



These results contradict the claim in [5] that the isocurvature perturbation is converted 
into an adiabatic perturbation after the equivalence. Our results are in agreement with other 
works (see [10] and references therein). To summarize, for the primordial isocurvature modes 
that remain outside the Hubble radius, Sk remains constant during the whole evolution, 
whereas the density perturbations 6k^ initially zero, grows like a"^ during the radiation era, 
then like a during the matter era, while remaining small with respect to Sk- Once the 
perturbations enter the Hubble radius both Sk and 6k grow rapidly (see Fig. 1). 

The same analysis for a primordial adiabatic perturbation reveals that, in the initial era, 
6k grows like while Sk^ initially zero, evolves like 




(20) 
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During the matter era, 6k grows like a. We refer the reader to [10] for a more detailed 
discussion on adiabatic perturbations. 



III. THE COSMIC MICROWAVE BACKGROUND RADIATION 

We turn now to the relation between the perturbations, discussed in section 2, and the 
observed CMBR anisotropy. The observed CMBR is the image of the last scattering surface 
which arrives today to our eyes (or in fact to our radio antennas). The fluctuations in the 
observed temperature are the consequence of both the perturbations in the matter content 
of the universe at the time of the last scattering and the perturbations of the geometry in 
the regions between the last scattering and us. 

The temperature fluctuations due to the matter perturbations at the time of last scat- 
tering are intrinsic fluctuations. Since the radiation energy density is proportional to T"^ one 
always has: 

AT 1 Spr 



For an adiabatic perturbation [S = 0), (11) implies 



(23) 



6r = P-ra. (24) 

The last scattering occurs during the matter dominated era, i.e. when p^ < Pm^ hence 
~ pSra and therefore 

Ar\ 1 . . 

I ~ — (5, adiabatic perturbation, matter era. (25) 

T J int 3 

For an isocurvature perturbation, Sp^ = —Sp^- Hence during the matter era, 6^ can be 
neglected with respect to 6^ and (11) yields S ~ — (3/4)(5r. Therefore 

Ar\ 1 

— - j ~ S, isocurvature perturbation, matter era. (26) 
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The second source of the CMBR temperature fluctuations is the so-called Sachs- Wolfe 
effect [14]. It corresponds to the influence of the metric perturbations on the light rays during 
their travel between the last scattering surface and "our eyes". To express the Sachs- Wolfe 
effect we introduce the gravitational held t/j and the peculiar velocity v. The gravitational 
held, is given by the relativistic generalization of the Newtonian Poisson equation with a 
FRW background, 

^H^S = D\l^. (27) 

During the matter dominated era, as stated in section 2, the dominant solution for the 
density perturbation grows like (5 ~ a and therefore the gravitational held is constant in 
time. Ignoring the decaying mode, the peculiar velocity held is deflned as (see e.g. [10]) 

^=-t3tjj. (28) 

Using t/j and v one can write the Sachs- Wolfe contribution to the CMBR anisotropy as 
/Ar\ ^1 

-7^ (e) = - [V'em - V'o] + e. [Vem " ^^o] , (29) 

\ J- y SW o 

where e is the unit vector corresponding to the direction of observation on the celestial 
sphere. The subscript em means that the corresponding quantity is evaluated at the point 
on the last scattering surface that is observed today in the direction e. The subscript 
refers to the observer today. 

The combination of equations (25) and/or (26) with equation (29) gives us the total 
CMBR anisotropy observed today (at least for large angular scales) for any conflguration of 
the adiabatic and isocurvature perturbations, which is here specifled in terms of the functions 
t/j[r) and S{r) at the time of last scattering. 

The observed CMBR temperature fluctuations are generally decomposed in spherical 
harmonics : 

— (^,<^) = E E (30) 

/=1 ra = — l 
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where 6 (varying between and tt) and (f> (varying between and 27r) are the usual angle co- 
ordinates on the two-sphere. One can express the individual coefficients of the decomposition 
as 

«/™ = j d<f>,\nede^{e,<f>)Y;^{e,<f>), (31) 



In the particular spherically symmetric case, which we examine in the following, the 
temperature fluctuations depend only on the angle 6, and all the coefficients with m 7^ 
vanish ( [15]). The non vanishing dipole and quadrupole coefficients are 




where u = cos 6. 

In the expression (31) the temperature anisotropy can also be seen as a function of 
X and can thus be decomposed in terms of its Fourier modes defined according to (18). We 
obtain 




where 0^; is the solid angle corresponding to x. x is the norm of x and represents the 
comoving distance between the observer and the last scattering surface. To a very good 
approximation, x ~ 2[Hoao)~^ (the exact expression is given in (56)). Using then the 
identities 

^rk..- = iTrJ2i'ji{kx)Yun,)Y;jnk) (36) 

l,m 
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and 

f J011,„,(0)l^:„^(0) = S^.iJ^,^,, (37) 

one finds 

= 4^/ J^2''Mkx)Y^:a{^,) {^)^- (38) 

For kx > /(/+ l)/2, 

~ 7- sm{kx - /7r/2). (39) 

In the expression (38) tlie small scales are suppressed due to tlie presence of tlie factor 
ji[kx). A rougli estimate of tlie coarse graining scale is x for the dipole and x/3 for the 
quadrupole. In any case the modes that contribute to (32) and (33) are large scale modes, 
whose wavelength is far larger than the Hubble radius at the time of last scattering. The 
latter corresponds to an angular scale of 1° today. In practice we are thus allowed to work 
with the perturbations smoothed on a scale of the order of a few times the (comoving) 
Hubble radius at the last scattering. On one hand, these smoothed perturbations will not 
change the results in (32) and (33) at the notable exception of the dipole term due to —e.VQ. 
On the other hand, these smoothed perturbations contain only Fourier modes that remain 
outside the Hubble radius until the last scattering and are thus far easier to handle, as shown 
in the previous section. 

To conclude this section we recall here that the satellite COBF [16] has measured the 
dipole and quadrupole components of the CMBR anisotropy to be: 



D 



Q 



1 1 

A T- E ^2x10-3, (40) 

\ m = — 1 



1 2 

J2 |a2™P^5xlO-^ (41) 



ITT ^_ _2 



\|4 

In fact the directly measured dipole is of the order D ~ 1.23 X 10""^ but we have extrapolated 
this result, as is usually done, to get a dipole measured with respect to the Local Group 
frame [17]. 
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IV. THE TOLMAN-BONDI MODEL 



The purpose of this section is to recover the numerical results of PP by an analytical 
method which is based on the linear theory presented in the previous sections. Indeed, 
although PP treated the complete non linear Tolman-Bondi (TB) problem, one can see from 
the order of magnitude of their quantities that one can study the same problem within the 
linear approximation. We begin, therefore, with the complete TB solution and we linearize 
it around the background solution of a flat FRW model. 

The TB solution is given by the metric 

ds"^ = -dt"^ + X^ir, t)dr'^ + R^{r, t)dn\ (42) 

where the coordinate r is a comoving coordinate attached to an element of the dust fluid. 
The solution is characterized by two arbitrary functions of r. The first, W{r)^ expresses how 
bound is a given shell r with the binding energy being W{r) — 1. The second ts{r) describes 
the time in which the radius of a given shell vanishes and, following PP, we will consider 
here only solutions with ts = 0. We define the gravitational mass as 

mb{r) = r 3r'^W(r)dr. (43) 
Jo 

The metric functions X and R satisfy then 

ti' = H,-'(r) - 1 + (45) 

R 

where a dot denotes a partial derivative with respect to the time coordinate t and a prime 
denotes a partial derivative with respect to the radial coordinate r. The mass density /?„(?", t) 
is given by the expression 
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Our background model is a dust dominated flat FRW universe, corresponding to the 
particular TB solution with W{r) = 1. One can easily solve (44) and (45) in this case, and 
with the initial condition i?(r, 0) = 0, one flnds 

X = a(t), R = a{t)r. (47) 

The FRW scale factor a[t) is explicitly given by 

a{t) = afl\ a = {^I2fl\ (48) 

Consider now a small deviation from this flat FRW model: 

W{r) = 1 + 8W{r). (49) 

Then the linearization of equations (44) and (45) yields 

8X = 6R' - aSW, (50) 

and 

hrSR + ^SR = SW+— [ Sr^SWdr. (51) 
a/ ar J 

The second term on the right hand of this last equation decays like with respect to the 
flrst one. Dropping it, one can easily flnd the dominant solution: 

,R = ^^E;./.. ,.,2) 

10a r ^ ^ 

The dominant solution for 8X follows immediately: 

,x=^f^Ey<./3. ,53) 

lOa V r ; ^ ^ 

One can also obtain the dominant solution for the matter energy density (5^: 

,„._^(!^2l,/3._^,i,.,-.ME):. ,54, 
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Substitution of expression (54) for 6^ in the Poisson equation (27) and using (48) yields 
an expression for the gravitational potential: 

V'(r) = -— [ dr ^''^^^ . (55) 
5r J r 

t/j is time independent as expected. 

It remains to calculate the temperature anisotropy for this spherical symmetric model. 
Consider an observer located at Tq at time to. There is a preferred axis which links this 
observer to the center r = 0. Since there is an axial symmetry around this axis, one can 
restrict oneself to the meridional plan. One then denotes 6 the angle between the preferred 
axis and the direction of observation [6 varies between and tt). Let d be the comoving 
distance between the observer and the last scattering surface: 

d = 2H-\-^ [l - (1 + ^)"'/'] = ^a-Hl^^ [l - (1 + ^)"'/'] , (56) 

where z is the redshift corresponding to the last scattering surface (for the numerical appli- 
cation, we shall take z = 1000). Then the radial distance of a point of the last scattering 
surface corresponding to the angle of observation 6 is given by 

= d^ sin^ 9 + {ro + dcos9 f . (57) 

Finally, using (29) and (28) the temperature anisotropy of the CMBR is given by 

(58) 

If the function AT/T varies sufficiently slowly one can expand it in a Taylor series and 
evaluate analytically the integrals (32) and (33). Starting with the Sachs- Wolfe part of the 
temperature fluctuation, one can see that there is a cancellation in the non symmetric terms 
(i.e. the dipole terms) at the order {d^ jr^) and that one must therefore go to the following 
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order [d /r^)^ which also implies that terms proportional to u will appear. One finally 



obtains 



— (u) ~ Diswu + D2SWU + Qswu 
1 J sw 



(59) 



wi 



ith 



sw 



1 di^ rld^tP 



(60) 



D 



2SW 



1 dtp r^d^tj) r-Qd^tp 



2 dr^ 



(^0) + 



6 dr^ 



J2 (d 



^-^r^ + ^a-H]l\l + z)-^r^]+Or-A (61) 



1/2 



J4 



and 



QsFP' — d 



d^ijj^ 
dr 



dr 



■^ + a-'4"(l + .)-"A+og' 



(62) 



The intrinsic temperature fiuctuation is simpler to evaluate. If we note [AT /T)int = /, 
where / is either 6/3 in the adiabatic case and —S/3 in the isocurvature case, then the 
Taylor expansion yields 

'AT' 



T 



int 



[u) ~ Diu + Qiu 



(63) 



wi 



ith 



Di = —{rQ)d + 0\^ 



dr 



2 ; ' 



(64) 



*4 



^0 - To ^ ro 



.o(C 



(65) 



Finally the total dipole and quadrupole, expresses in terms of the above coefficients, are 



1 / 3 
T) = -j^[Di^ Disw + -D2SW 



(66) 



anc 
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Q = —^{Q, + Qsw) . (67) 

The expression seems similar lor both the isocurvature and adiabatic perturbations. 
However, the Poisson equation (27) with the expression (56) shows that 8 is always ol the 
order (J/ro)^(l -\- z)~^ times Therefore the intrinsic dipole contribution due to 8 is smaller 
than the Sachs- Wolfe dipole by a factor of the order (1 -\- z)~^ and is thus always negligible. 
But when there is an isocurvature perturbation, one sees, in contrast, that the intrinsic 
dipole will be in general the dominant term. 



A. PP results 



We are now in position to compute analytically the results of PP. Their specific model 
corresponds here to 

1 1 - r2 , 



8W{r) 



and 



S{r) 



2 1 + r 



1 



2 ' 



2 1 + r2 

It turns out that this expression for 8W{r) yields an explicit expression for ?/'(r): 

3 
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1 — r ln(l + r^) 

r 1 

1 + 2 r 



The numerical results are the following (with to = 10 



(68) 



(69) 



(70) 



Di{iso) ~ -1.47 X 10-^ Disw ^ -1-26 x 10-^ D2SW ^ 5.09 x 10" 



(71) 



Qiiiso) ~ 2.58 X 10-^ Qsw ^ -6.94 x 10'^. (72) 

For the model without isocurvature perturbation, one finds D ~ 1.0 X 10"'' and Q ~ 
2.1 X 10~^, whereas for the model with isocurvature perturbation, one finds D ~ 8.5 X lO""* 

17 



and Q — 1.9 X 10~^. These numerical results should be compared with the results obtained 
by numerical integration of the light rays, and given in the figures 3 and 4 of PP. Our results 
here are limited to Tq = 1 since we have assumed a fiat space from the beginning. The 
numerical values given above correspond to to = 10~^ but the corresponding values for a 
different to can be obtained immediately by noticing that the dependence on to of Dj is 

1/3 

due to an overall multiplicative term to whereas Disw and D2SW are proportional to to, 

2 /3 

and Qsw and Qi proportional to to . Comparison with PP shows a good agreement, thus 
confirming the conclusions of PP, although there are small discrepancies between the precise 
numerical values, which we cannot explain. 

One could have argued against the results of PP that their two perturbations are a priori 
completely independent. The question arises what will happen if one considers an initial set 
of primordial perturbations and let it evolve in time until the time of last scattering. Will 
it be possible to reproduce similar results with these more stringent conditions? In [5], it is 
argued that it will be impossible to recover PP results. What we show in the next subsection 
is that indeed we can recover the same behavior. 



We turn now to consider the extreme case where primordial isocurvature perturbations 
on extremely large scale (i.e. scales much larger than the horizon) are the only source of the 
observed dipole and quadrupole. We denote by S^. the modes of the primordial isocurvature 
perturbation [k < aoHo). It follows from the analysis of Section 2 that: 



where Is stands for the last scattering. The primordial isocurvature perturbation has also 
generated a energy density perturbation given by 



B. Primordial isocurvature perturbation 




(73) 




(74) 
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and a corresponding gravitational potential: 



(75) 



We choose, as an example, a primordial isocurvature perturbation of the form 





(76) 



1 + [r/rsY 



One can adjust the two parameters of the perturbation, namely the amplitude and the 
wavelength r^, to reproduce exactly the measured dipole and quadrupole. One finds ~ 
2.94 and ~ 2.80, which means that the wavelength of the perturbation must be roughly 
150 times larger than the Hubble radius today. Note that such adjustment is possible only 
for isocurvature modes. This would be impossible with an adiabatic mode, as will be shown 
clearly in the next section. 



Before turning to observational implications, we wish to extract, in this section, the 
essential arguments which explain why the isocurvature perturbation and the adiabatic per- 
turbation produce such different results. To show this we consider a plane wave perturbation 
and obtain rough estimates for all terms involved. We show in particular that the ratio D/Q 
for a very large scale perturbation is inverted when one goes from an adiabatic to an isocur- 
vature perturbation. 



We have already emphasized that the main contribution to the dipole or even to the 
quadrupole arises from large scales. We shall restrict, therefore, our analysis to the Fourier 
modes outside the Hubble radius (at the time of the last scattering). The Fourier transform 
of the relativistic Poisson equation (27) reads 



V. INTERPRETATION OF THE RESULTS 



A. Adiabatic perturbations 
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Assuming a monochromatic perturbation of the form 8 ~ (^^e'*^^, the intrinsic dipole contri- 
bution is, therefore, of the order 

The dipole contribution of r/j/S is given by ( [11]) 

D{^-i^] = -<v{to)-v{tu)>, (79) 

where the average is taken over the comoving volume defined as the intersection of our past 
light-cone with the last scattering hypersurface 

<v{t)>=Vi;^ [ d^xv{t,x). (80) 

J Is 

Finally, we recall that the peculiar velocity evolves like t/a. Thus, the peculiar velocity 
at the last scattering is negligible with respect to the peculiar velocity today and the main 
contribution to the Sachs- Wolfe dipole for adiabatic perturbations is 

Dsw ^ e. [vo- <v> (to)) . (81) 

This result corresponds to the standard statement that the dipole is due to the relative 
motion of the Earth. Note that this calculation gives a precise definition of the relative 
velocity of the Earth and in particular with respect to which frame. 
Using a Taylor expansion within the integral in (80) one obtains: 

Dsw ^ e. [vo- <v> (to)) ~ (-YTHI "fc(^o) ~ ( — TT ) (^2) 

Vq corresponds, in this formula, to our peculiar velocity induced only by the single very large 
scale mode under consideration. If one takes into account the contribution of small scale 
modes to our peculiar velocity then one sees that the net measured velocity (and hence the 
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measured dipole) will be dominated by the contribution of the sub-horizon modes. Moreover, 
comparison of (78) with (81) and (82) shows that Di{6} ~ (1 + z)~^Dsw and therefore the 
total dipole is D ~ Dsw- 

The expected dipole can be evaluated given a power spectrum (the adiabatic pertur- 
bations are supposed to be distributed like a Gaussian random field). Using the power 
spectrum of t/j defined by 

(V-.^Vi,) = 27r' k-'V4k)6{k - k'), (83) 

we find, taking into account only the (dominant) term e.Uo, 

167r f dk f P \ _ , , , 

For the most common spectrum, namely the scale invariant Harrison- Zeldovich spectrum, 
is a constant and we express it as a function of the expected quadrupole 
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= -E^, (85) 

TT 

with E2 = (|a2mP)- The contribution to the rms-dipole of the fiuctuations of scales larger 
than a given scale is obtained by integrating the above expression (84) with a lower cut-off 
k. Hence 

("V^ ^ ^i^r^- - S.^. (86) 

This gives for the peculiar velocity 

v.rms ^ 250A;m s'^ (Qr^s-psin = l)/17i^K) (A/50/i-^ Mpc) , (87) 

where h is the Hubble constant in units 100 km s~^Mpc~^ and A is the smoothing scale 
for the fiuctuations (the relation between Qrms-ps{n- = 1) such as it is defined in [16] and 



E2 is E2 = y 4:Tr /5[Qr'ms-PS /Tq)^ where To = 2.73A^ is the CMBR monopole temperature). 
Instead of using a cut-off in the Fourier space, one may wish to introduce a cut-off in the 
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real space, i.e. use a top hat window function, in which case the expression (86) should be 
multiplied by the factor V 4.5 ~ 2. This relation shows that our velocity measured with 
respect to a given frame is inversely proportional to the distance of this frame from us. 
The CMBR quadrupole is given by the quadrupole part of r/j/S which is of the order 

The expected ratio between the dipole and quadrupole contributions arising from a single 
mode very large scale adiabatic perturbation is, therefore, of the order 

This is less than unity, by definition, hence the observed dipole and quadrupole cannot be 
explained by such a perturbation. In the context of adiabatic perturbations, the origin of 
the dipole must be only peculiar velocity. 

Note finally that, since aH = (1 + zY^^aoHo and aH ~ tg^^^ in the matter dominated 
era, expressions (82) and (88) show that D ~ to, Q ~ tg^^ and {D/Q)ad ~ tg^^. This is in 
agreement with the numerical dependence observed by PP. 



B. Isocurvature perturbations 

We turn now to isocurvature perturbations. Whereas the dependence on the gravitational 
potential t/j of the Sachs- Wolfe term due to a primordial isocurvature perturbation is the same 
as that due to a primordial adiabatic perturbation, the intrinsic contribution is drastically 
different. In particular the dipole of the intrinsic anisotropy, 

is dominant with respect to the Sachs- Wolfe dipole for scales larger than the Hubble radius 
today. The total quadrupole, on the other hand, is comparable to the adiabatic one: 
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Q^^o ~ (-^772) V'fc- (91) 
\"o-"o / 

Combining the last two expressions, we find that the expected ratio between the dipole and 
quadrupole is 

( D\ ( k 



(92) 



\Q J ■ \ (IqHq I 

\ ^ / ISO \ / 

Once more, a comparison with the numerical behaviours observed in PP is instructive. The 

2 /3 

quadrupole Q has the same form as in the adiabatic case: Q ^ . The dipole is difi'erent 

1/3 

and it follows from the above expression that it evolve like D ~ . The ratio behaves now 

— 1/3 

like {D IQ)iso ~ . All these dependences are confirmed by the numerical results of PP. 

C. Continuous isocurvature spectrum 

So far we have dealt only with ultra large scale monochromatic perturbations. One can 
wonder what will be modified when one considers a continuous spectrum instead of a single 
mode. We examine here only a pure isocurvature primordial spectrum. The case of an 
adiabatic power spectrum is extensively treated in the literature (see e.g. [10]). 

We consider primordial isocurvature perturbations that are described by a homogeneous 
and isotropic Gaussian random field, which is completely specified by its power spectrum 
Vs{h): 

{Sj:Sl,) = 27r'k-'Vs{k)S{k - k'). (93) 

In particular one can compute the variance of the distribution of the multipoles as a function 
of the power spectrum Vs'- 

= {\aim?) = yJ frs{k)jf{2k/aoHo), (94) 

where it is assumed that the intrinsic contribution is the dominant one in the temperature 
anisotropy. The dipole and quadrupole correspond roughly to Ei and E2 respectively. We 
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now assume that the power spectrum is a power law Vs{k) ~ Ak^^ which is a standard 
assumption made in cosmology, and moreover that this power spectrum as an upper cut-ofF 
k. 

If the cut-ofF k is such that k <^ Oo-f^o, then it is legitimate to use the approximation 

x' 

One can then calculate explicitly the dipole and the quadrupole: 

D ~ U-^+'y\aoHo)-\ Q ~ U-^+'y'{aoHo)-\ (96) 

The ratio between the dipole and quadrupole due to a very large scale (power-law) power 
spectrum should thus be of the order 

D ( k 



(97) 



Q \aoHo J 

which is of the same order as in the case of a monochromatic perturbation. 

If there is no cut-off in the power spectrum or if the cut-off is smaller than the Hubble 
radius today, i.e. k > ao-ffo, then the dominant contribution in the integral (94) both for 
1=1 and 1 = 2 will come from the wavelengths roughly of the same order than the Hubble 
radius (today), as explained at the end of section 3. Therefore the corresponding dipole and 
quadrupole should be of the same order of magnitude in this case. Beware that the dipole and 
quadrupole here are computed by taking into account only the intrinsic contribution. For 
the scales smaller than the Hubble radius (today), the primordial isocurvature modes have 
produced adiabatic perturbations with corresponding gravitational potential and peculiar 
velocity field. The dipole can therefore be dominated by the Doppler effect due to our 
peculiar velocity, as in the standard interpretation of the dipole, and we thus find that the 
dipole can be large relative to all other multipoles, even in this case. 

Finally we must mention the hybrid possibility that the observed dipole could result 
of a combination of a Doppler effect (due to either adiabatic or isocurvature primordial 
perturbations) and of an intrinsic ultra large scale isocurvature contribution. 
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VI. CONCLUSION 



We have shown that it is possible that the observed CMBR dipole, or a significant frac- 
tion of it, has a non Doppler origin. This can happen if there is a suitable ultra large scale 
(typically beyond 100 times the size of the present Hubble radius) isocurvature perturbation. 
That is an ultra large scale fiuctuation in the ratio of photons to baryons. This fiuctuation 
will induce predominantly a dipole component in the observed horizon, without inducing 
higher order multipoles. This fiuctuation should not be accompanied by isocurvature per- 
turbations on smaller scales (between IOOHq^ and O.IHq^), because these would induce 
higher order moments which would then be comparable to the dipole. 

It is clear from the above analysis that if the observed CMBR has a non Doppler origin 
then there should be a unique mechanism that would produce these very large scale isocur- 
vature perturbations and will distinguish them from the rest of the power spectrum. A 
priori there are several possible origins for these ultra large scale perturbations. They could 
have been produced during an infiation era: the model of infiation (see e.g. [18]) must then 
include several scalar fields in order to allow for isocurvature perturbations in addition to 
the always present adiabatic ones. Another explanation would be that these perturbations 
are the remnants of the preinfiationary epoch of the universe, as was suggested by Turner 
[4]. Indeed, if the duration of infiation is slightly more than what is needed to solve the 
"horizon" problem, then the scales that were larger than the Hubble radius at the onset of 
infiation would be today also larger than the Hubble radius but not by a large amount. 

Is it possible to distinguish between an isocurvature dipole and a dipole due to our 
peculiar velocity? Luckily, there is a clear direct observational test. We have seen in the 
last section that the dominant contribution to our peculiar velocity arises from small scale 
modes and hence it should converge to the same velocity when it is measured relative to 
different distant frames. Thus a peculiar velocity dipole will induce the same dipolar pattern 
in other background fields, such as the X-ray background or 7-ray bursts which are located 
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at 2; 1 — 2. Depending of the power spectrum even nearer frames such as optical galaxies, 
IRAS galaxies, distant supernovae and Abell clusters should display similar peculiar motion 
pattern. Note however, that for distances smaller than the horizon the intrinsic contribution 
to the dipole might not be negligible and fluctuations in the density of sources should be 
included [19]. A convergence of all those peculiar velocities will support the Doppler origin 
of the CMBR. A final confirmation should arise if the observed peculiar motion is consistent 
with the expected r.m.s. value of this quantity, given by Equation (86), as calculated from 
the observed power spectrum of the matter fiuctuations (one has of course to be careful here 
about cosmic variance). If this interpretation is confirmed then the measured quadrupole 
shows that the Universe is homogeneous at least on scales that are larger by 10^ than the 
current horizon. This will immediately rule out the existence of significant ultra large scale 
isocurvature perturbations and cosmological scenarios that produce them. 

If, on the other hand, the observed dipole is due to extremely large scale isocurvature 
fiuctuations, it should not have any corresponding signature on small scales. We expect, in 
this case that the observed dipole relative to the nearer frames, mentioned earlier, will still 
converge, but now to a different velocity in both magnitude and direction than the velocity 
implied by the CMBR dipole. We should point out that failure of the peculiar velocity to 
converge on those nearer scale would imply that the primordial power spectrum has some 
peculiar behavior on intermediate scales, a behavior that causes the integral in equation (94) 
to fiuctuate. 

The observed CMBR dipole implies a velocity of the local group of 627 it 22 km s~^ and 
it points towards the galactic coordinates (/ = 276° it 3°,6 = 33° it 3.°). The magnitude 
of this velocity is larger than the expected r.m.s. value for a Harrison-Zel'dovich spectrum 
normalized by COBE, which is several hundred km s~^. The measured dipole in the X- 
ray background [20], which arises from sources at 2; 1.5 is within the statistical errors. 
The dipole has not been measured relative to any other sources at comparable distances. 
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However, it has been measured relative to nearby galaxies [21], IRAS galaxies [22] and distant 
supernovae [23] which are all at 2; < 0.03. The COBE observations are within the statistical 
errors of all those measurements. This suggests that we are observing the convergence of the 
dipole on smaller scales. The observed dipole in distant Abell clusters [24] whose magnitude 
is 561 ± 284 km s"^ towards (/ = 220° ± 27°, b = -28° ± 27°) is, however, inconsistent with 
the CMBR dipole. It is also inconsistent with the distant supernovae dipole [23] which is 
measured relative to objects at the same distances. Hence we can conclude that at present 
the observational situation tends towards the conventional Doppler origin but the situation 
is inconclusive yet. Further measurement in the future of dipole relative to additional frames 
or refinement of current measurements should provide a conclusive answer in the future. 
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FIGURES 



Fig. 1: The amplitude of the entropy perturbation (Sk) and the density perturbation 
(8k) as a function of \og{a j apq) , for an initial pure isocurvature perturbation with an initial 
amplitude 10~^. Three different wavelengths are shown. The first (solid curve for S and 
dotted curve for 8) remains always larger than the horizon. The second (short-dashed curve 
for S and long-dashed curve for 8) enters the horizon at a ^ l.baeq, which is marked in the 
figure by a square. The third perturbation (dashed-dotted curve for S and curve da.shed- dotted 
line for 8) enters the horizon at a ^ 0.05aeg, which is marked by a triangle. 

Fig. 2: The quantity [aH j kY 8k as a function oflog[a/aeq) for a primordial pure isocur- 
vature perturbation with five different wavelengths. The two long wavelengths (dotted curve 
and long-dashed-short-dashed curve) overlap. These modes do not enter the horizon and they 
approach the asymptotic limit (straight solid curve) 2/15. An intermediate wavelength per- 
turbation (.short dashed curve) deeps slightly before going up after it has entered the horizon 
at a ^ apq. The same behaviour is seen in the two .short wavelength perturbations (solid curve 
and long-dashed curve) which deep first and then grow rapidly after entering the horizon at 
a < O.laeg. 
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